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ABSTRACT 


This  report  contains  the  results  of  the  research  project  entitled: 
"Microcracking  and  Toughness  of  Ceramic -Fiber/Ceramic -Matrix  Composites 
Under  High  Temperature"  supported  by  the  AFOSR.  Microcracking  mechanisms 
and  toughness  of  Nicalon  (SiC)/SiC  composite  at  elevated  temperature  are 
studied  analytically  and  experimentally.  First  the  fiber  distribution 

patterns  in  the  ceramic  composite  are  determined  by  observing  the  specimens 
under  optical  and  scanning  electron  microscopes.  Thus  the  effect  of  fibers 
and  fiber  interactions  on  the  microcrack  propagation  are  investigated 
analytically  through  the  single-fiber,  the  two-fiber,  and  the  ring  models. 
Monolithic  SiC  specimens  are  tested  under  varying  temperature  to  determine 
the  effect  of  temperature  on  the  toughness  of  the  matrix  material.  The 
Nicalon/SiC  composite  specimens  are  then  tested  at  various  temperatures. 
The  combined  effect  of  temperature  and  fibers  on  the  toughness  of  the  matrix 
is  expressed  by  introducing  the  concept  of  "apparent  fracture  toughness". 
The  experimental  results  indicate  that  for  the  Nicalon  /SiC  composite  the 
"apparent  fracture  toughness"  decreases  with  local  volume  fraction  of  fibers 
and  temperature.  An  analytical  model  to  predict  this  behavior  is 
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1. 


Introduction 


As  is  well  known  ceramic  matrix  composites  have  recently  attracted  a 
great  deal  of  attention  due  to  the  promise  they  hold  for  developing  high- 
temperature  resistant  materials  with  high  toughness.  Up  to  now  the 
structural  applications  of  ceramics  and  ceramic  composites  have  been  limited 
becaqse  of  the  fact  that  their  brittleness  makes  them  an  unlikely  candidate 
in  applications  where  high  toughness  is  required.  But  the  combination  of 
advanced  manufacturing  techniques  and  new  ceramic  composites  appears  to 
satisfy  both  conditions  of  resistance  to  high- temperature  oxidation  and  high 
toughness.  The  current  study  focuses  on  one  such  material,  namely, 
Nicalon/SiC  composite  manufactured  by  Amercom,  Inc.,  Chatsworth,  CA. 

It  has  long  been  recognized  that  the  low  toughness  of  ceramics  results 
from  the  presence  of  microcracks  or  voids  In  the  material.  These 
microcracks  may  develop  due  to  inclusions,  thermal  expansion  mismatch,  phase 
transformation,  or  thermal  expansion  anisotropy.  Whatever  the  initial 
agent,  under  sustained  mechanical  and/or  thermal  loading,  these  microcracks 
may  coalesce  into  a  critical  flaw,  causing  the  failure  of  the  material.  In 
fiber  reinforced  ceramic  composites  the  failure  mechanism  is  much  more 
complicated  because  of  the  presence  of  fibers.  For  example,  under  uniaxial 
tension  fiber  matrix  microcracks  may  appear.  At  sustained  loading  these 
cracks  may  join  other  cracks  to  form  a  dominant  flaw  or  may  be  deflected  at 
the  fiber/matrix  interface.  It  appears  that  the  fiber/matrix  interface 
plays  an  important  role  in  the  overall  toughness  and  strength  of  the 
composite.  For  instance,  a  relatively  weak  interface  may  help  deflect 
matrix  cracks  while  a  strong  interface  may  lead  to  fiber  fracture  and  the 
eventual  catastrophic  failure  of  the  composite.  The  literature  on  factors 
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affecting  the  toughness  of  ceramics  and  ceramic  composites  is  extensive. 
The  microstructural  dependence  of  fracture  parameters  of  ceramics  are  given 
in  [1-12].  Particle  and  grain  size  effects  are  studied  in  [13-17]. 
Microcracking  in  particulate  ceramic  composites  due  to  thermomechanical 
stress  is  investigated  in  [18]. 

In  this  work  microcracking  in  Nicalon/SiC  ceramic  composite  is  studied 
for  thermomechanical  loading.  Monolithic  SiC  (matrix  material)  and 
composite  specimens  are  first  observed  under  microscope,  and  the  fiber 
distribution  patterns  are  identified.  Thus  analytical  studies  relevant  to 
typical  fiber  distribution  patterns  (e.g.,  a  single  fiber  embedded  in  an 
infinite  medium,  two  isolated  fibers,  a  ring  distribution  of  fibers,  etc.) 
are  performed  and  the  consequences  on  matrix  microcracking  are  discussed. 
Next  the  toughness  of  matrix  material  (SiC)  is  determined  under  various 
temperatures  (up  to  800  C)  using  the  micro- indentation  technique  [19,20]. 
The  effect  of  fibers  on  the  "apparent  fracture  toughness"  of  matrix  is 
determined  experimentally  again  using  the  micro- indentation  technique.  A 
theoretical  model  predicting  the  microcracking  behavior  is  then  developed. 
Finally,  the  debonding  strength  of  the  interface  at  the  initiation  of 
debonding  is  determined  by  a  method  combining  experimental  results  and 
finite  element  calculations. 

2 .  The  Material:  Nicalon/SiC  Composite 

The  material  studied  is  Nicalon/SiC  purchased  from  Amercom,  Inc. , 
Chatsworth,  CA.  The  properties  of  the  const.' tuents  are  given  in  the 
following  table: 


U 


Young's  Modulus 

Coefficient  of  Thermal 

E  (nsi) 

Exoansion  *  (in/in/  C) 

Fibers 

(Nicalon) 

26xl06 

3. 1x10' 6 

Matrix 

(SiC) 

55xl06 

4. 3x10* 6 

It  should  be  noted  that  this  material  is  a  reversed  composite.  That  is, 
stiffness  of  the  matrix  (SiC)  is  greater  than  that  of  the  fibers  (Nicalon) . 
The  volume  fraction  of  fibers  is  approximately  0.35  and  the  void  volume 

fraction  is  about  0.15.  The  composite  is  manufactured  using  the 

chemical  vapor  deposition  (CVD)  method. 

Using  the  information  obtained  from  preliminary  testing  and  considering 
the  size  of  the  heating  chamber  of  the  furnace,  the  test  specimens  were 
designed  as  follows: 

12  pc.  of  (1/2) "x(l/2) "x(l/4) "  3-D  braided  Nicalon/SiC  specimens 

infiltrated  by  CVD  to  a  minimum  density  of  2.3  g/c.c. 

6  pc.  of  (1/2) ”x(l/2) "x(l/4) "  silicon  carbide  coated  graphite 
specimens  coated  on  one  face  with  SiC  by  CVD  to  a  minimum 
thickness  of  0.020  inches. 

The  first  set  of  specimens  is  used  to  determine  the  effect  of  fibers  on 
matrix  microcracking  and  the  shear  strength  of  the  fiber/matrix  interface. 
The  second  set  is  used  to  determine  the  fracture  toughness  of  the  matrix  at 
varying  temperatures. 

An  examination  of  the  composite  specimens  under  microscope  reveals  that 
there  are  basically  four  fiber  distribution  patterns:  (a)  an  isolated  fiber 
(Photo  1),  (b)  two  isolated  fibers  (Photo  2),  (c)  a  circular  array  (ring)  of 
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fibers  (Photo  3),  and  (d)  closely  packed,  randomly  distributed  fibers.  The 
microcracks  could  be  in  the  matrix  or  at  the  interface.  Thus  the  analytical 
studies  presented  in  the  next  section  are  based  on  models  simulating  the 
observed  configurations. 

3 .  The  Theoretical  Studies 

The  analysis  for  each  specific  model  considered  is  outlined  below. 

(a)  The  single-fiber  model 

This  model  is  designed  to  simulate  situations  seen  in  Photo  1.  It  is 
assumed  that  a  single  fiber  is  embedded  in  an  infinite  matrix.  The 
microcracking  configurations  studied  are  shown  in  Fig.  1. 


(Muff/ X 

matrix 

f rr) 

■<p^ 

-/■be  r 

i  rt  i  °v 

! 

(a)  Radial  Crack  in  Matrix  (b)  Interface  Crack 


FIGURE  1. 

First  the  stress  fields  for  the  uncracked  configurations  are  computed 
in  closed  form  for  temperature  rise  as  well  as  mechanical  loads.  Consider  a 
ceramic  fiber  of  radius  R  embedded  in  an  infinite  ceramic  matrix  as  shown  in 
Fig.  2.  If  the  temperature  is  increased  by  T,  the  stresses  in  the  fiber 
and  the  matrix  can  easily  be  found  by  considering  a  uniformly  loaded 
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cylinder  and  a  plane  with  a  hole  and  matching  the  displacement  at  their 
common  boundary.  After  some  algebra,  we  obtain: 


for 

r  <  R  - 

r 

arr 

—  a 

(in 

the 

fiber) ; 

(1) 

>  aee 

-  a 

for 

r  >  R  • 

«■ 

°xx 

■  °66 

-  cr  (— )  2 

\r/ 

/  R  \  2 

-  -a (  ) 2 

vr 

(in 

the 

matrix) ; 

(2) 

where 


a 


an/1+lO  ' 
m _ m _ t _ i 

A 


■AT; 


and 


A 


l+i/ 


m 


E 

m 


+ 


(3) 


(4) 


Here  (Em,  am)  and  (E^,  a are  the  Young's  moduli,  the  Poisson's 

ratios  and  the  coefficients  of  thermal  expansion  for  the  matrix  and  fibers, 
respectively.  Assuming  AT  >  0,  i.e.  heating,  and  considering  that 

one  may  conclude  immediately  that  if  am<  a^,  (a0 $ )matrix  *-s  tensile  and  may 

induce  radial  cracks  in  the  matrix.  On  the  other  hand  if  a  >  a^, 

m  f 

^rr^matrix  tenstle  with  the  possibility  of  debonding  along  the  fiber- 

matrix  interface.  In  practice  of  course,  in  addition  to  the  thermal 
stresses  described  above,  one  has  to  include  also  the  mechanical  stresses 
due  to  external  loading  and  the  residual  stresses  induced  in  the  material 
during  the  manufacturing  process. 
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Radial  crack  in  the  matrix.  Consider  the  radial  crack  geometry  shown 
in  Fig.  1(a).  The  formulation  of  this  problem  for  uniaxial  loading  is  given 
in  [21].  Here  instead  of  external  loads,  the  crack  surface  will  be  loaded 
with  the  residual  stresses  derived  in  the  previous  section  and  expressed  by 
Eq .  (2).  Referring  to  [21],  after  some  simplification  the  formulation  of 
the  problem  in  terms  of  the  crack  surface  displacement  derivative  f(t)  can 
be  written  as: 


c+R 

R 


fin 

t-x 


dt  + 


c+R 

[kiis(x, t)  +  kllf (x, t) ] f (t)dt 


*(3+1) 

m 


2m 


p(x) , 


R  <  x  <  c+R; 


m 


with  the  single  valuedness  condition: 


(5a) 


c+R 

f (t)dt  -  0;  (5b) 

R 


where 


.  A  rn  4s  s(s2-R2)  s3  (s2  -R2 ) 2  . 
Al  lu'  t  x2(t-s)2  '  R4t(t-s)3  J  ’ 


AlR2  2  1  3R2 


kHf(x,t)  =  x2  x"I_2 1  ~  tx 2 


C-K7T-73)  -  [M(E,+1)-1] 


1  -m  _ 

11  ’  1+mH  ’ 

m 


— c  -  nui  n  2 

f  m  R2 

A,  -  - - ;  s  -  — 

*  ^,^+m  x 


m  -  — ; 
“m 


M 


m(Em+1) 


(Ef+m) (Ef-l+2m)  ’ 


E  -  3-4i/  ; 
m  m 


~f  ”  3'4t/f  ’  Mm  =  2(l+u  )  ’  Mf 


m 


2(l+«/f) 


(6) 
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For  the  problem  undar  consideration,  the  crack  surface  will  be  loaded  with 
the  negative  of  residual  stress  a.,  given  in  Eq.  (2).  Thus, 


P(x)  =  0~.  (7) 

To  solve  the  singular  integral  equation  (5a),  the  following  normalization 
is  used: 


X  -  fp 

+  R  + 

c 

2’ 

for  R  <  x  <  c+R 

and 

- 1  <  p  <  +1 ; 

t  =  §r 

+  R  + 

c 

2’ 

for  R  <  t  <  c+R 

and 

-1  <  r  <  +1; 

f(t)  = 

g(0 ; 

kUs(x't) 

-  Kns(p,r); 

p(x)  - 

q (p) ; 

kuf  (x,  t) 

Then  Eqs . 

'  +1 
-1 


(5a, b)  become: 


^dr  +  f 

T-p  2 


+1 


[KUs(p,r) 


"(2+1) 

m 

2u 


q (p) . 


+  Kllf(p,r) ]g(r)dr 

-1  <  p  <  +1; 


and 

*  +1 

g(r)dr  -  0. 

-1 


(9a) 


(9b) 


The  normalized  crack  surface  displacement  g (r)  is  singular  at  r  -  ±1  and  may 
be  written  as  [ 22 ] : 


g  (O 


G(r ) 

yr 7d+r)^’ 


(10) 
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where  G(r)  is  a  bounded  function  and  -1  <  f)  <  0  is  given  by  [21]  as: 


Icosnfi  +  ( A i  +  A 2 )  -  4A1(/?+l)2  =  0. 


(ID 


At  the  crack  tip  embedded  in  the  matrix  the  stress  intensity  factor  can  be 
defined  as : 


k(o)  -  “R+c72lx-(R+c))‘,yy(*'0)| 


matrix 


(12) 


Using  Eqs .  (5a),  (8)  and  (10)  and  after  some  lengthy  algebra,  we  obtain: 


2/f 


k(c)  - 

“m 


(13) 


Noting  that  -  Jnk(c) ,  then  the  strain  energy  release  rate  is  found  to  be: 


K?d-^)  _ 

m  m 


(14) 


or 


GI  -  E 

m  J  c 


(15) 


C  \\  ( c  ^  — 

Defining  —  -  e  and  noting  that  —  -  <7k'(e),  Eq .  (15)  becomes 
K  — 


ys 


GI  “  E  (1-^)R£[k'(£)]V’  (16) 

ra 

The  radial  matrix  crack  will  propagate  when  the  strain  energy  release  rate 
reaches  the  critical  strain  energy  release  rate  of  the  matrix,  i.e.,  when 
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we  obtain  the 


GI  -  (GIc)matrix*  Equatin&  Eq.  (16)  with 
critical  fiber  size  as: 


(G  ) 

Ic'matrix’ 


R 


(Gt  )  _  .  E 

Ic  matrix  m  _ 

7r(l-i/2)e[k'  (e)]2a2' 


(17) 


However,  the  crack  size 
determine  only  a  minimum 
g(e)  -  €[k'(e)]2  is  maximum. 

Rmm  Ic  matrix  m 

c  “  7r(l-x/2 )  [g(< )  ] 

\  'Jmax 


(or  c)  is  not  known  a  priori, 
critical  fiber  size  which 
Thus , 


Therefore,  we  can 
will  occur  when 


(18) 


Formula  (18)  can  be  used  to  conservatively  determine  the  critical  fiber  size 
for  radial  matrix  crack  suppression. 

As  an  example,  the  fibers  are  assumed  to  be  SiC  and  the  matrix  a  glass - 
ceramic  (lithium  aluminosilicate  or  LAS)  with  the  following  properties  [23]: 


E  -  85  GPa; 

m 

Ef-  200  GPa; 

(Gc> interface 

(GT  )  .  - 

Ic  matrix 


v  - 

m 


-  20 


<KIc> 


0.2; 

a  -  9x10'6/°C; 

m 

-6  „ 

0.2; 

of-  4x10  D/  C. 

Pa-m ; 

AT  -  -1000°C; 

2 

matrix 

( 1  -i/2 )/E  -  45. 

m  '  m 

(K.  )  _  . 

Ic  matrix 

Pa-m . 


2  MPaVni; 


or 


Equations  (9a)  and  (9b)  are  solved  numerically  for  the  material  combination 
defined  above  by  using  the  collocation  technique  described  in  [24,25]. 

c 

Figure  3  shows  the  variation  of  g(<r)  with  ~.  It  is  seen  that  the  curve 
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passes  through  a  maximum  around  [ g( e ) ]  =*  0.075.  It  must  be  noted  that 

the  variation  of  g(e)  shown  in  Fig.  3  is  for  the  specific  material 
combination  used  in  the  example.  To  assess  the  effect  of  the  elastic 
properties  on  the  value  of  g(e),  results  are  obtained  for  various  ratios  of 
E^/Ef  ranging  from  0.2  to  5.  These  results  are  shown  in  Fig.  4. 

As  can  be  seen  from  Fig.  4,  surprisingly  the  Em/E^  ratio  has  little 

effect  on  g(«).  Therefore,  one  curve  can  represent  all  practical  material 
combinations,  including  the  Nicalon/SiC  being  studied  here.  Using  the  value 

[g(0]  “  0.075  and  substituting  into  Eq.  (18),  we  obtain  Rm*n=«  148pm, 

U13X  C 

which  is  larger  than  the  actual  fiber  size  used  in  such  composites. 

Interface  crack  along  the  fiber/matrix  boundary.  For  this  case,  the 
crack  geometry  is  shown  in  Fig.  1(b).  It  is  well  known  that  at  the 
interface  crack  tip  the  stresses  and  strains  have  an  oscillatory  singularity 
[26-30]  and  the  stress  intensity  factors  cannot  be  defined  by  using  their 
classical  definition.  Instead,  we  will  determine  the  strain  energy  release 
rate  at  the  crack  tip  which  is  not  affected  by  this  peculiar  behavior.  The 
problem  does  not  appear  to  have  been  solved  for  thermal  stresses.  Using  the 
solution  given  by  Toya  [31]  and  the  superposition  principle  we  present  here 
a  solution  for  the  residual  stress  loading  expressed  by  Eqs.  (1)  and  (2). 
For  the  plane  loading  shown  in  Fig.  5,  the  strain  energy  release  rate  at  the 
interface  crack  tip  is  given  as  [31]: 

G  -  ^kRA0  ( 1+4A^  )jrNNsinaexp  [  2A0  (7r-a)  ]  ;  (19) 

where 
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k  - 

Ao 

N  - 

G0 


H0 


d 


-1 


Po 

T^’  »  “  ^f+Sf%)/("m+Vf):  “  ^m(1+Hf)/(^m+Vf) 

,  1+H  1+E- 

S<— “  +  „);  *o  - 

^-1  1-k 

(c0— J^)  +  JYS(Noo-Ta3)exp[i(2^-a)+2A0(Q-7r)]  ;  c0  -  G0+iH0  ; 


[Jo^-TJ  ^■(cosa+2^osina)exP[2A0(7r-Q)  ]  } 

-J(l-k)  (1+4A§)  (Nm-Too)sin2acos2^]/ 

{ 2-k-k(cosa+2A0sina)exp[2A0  (*-a) ] } ; 
[J(l-k)  (1+4A§)  (Noo-T#o)sin2asin2<^ 

4^i  ie 

+  l+(cosa+2A0sina)exp[2A0  (n-a)  ]  }  ]/ 


m 


k{l+(cosQ+2A0sina)exp[2A0(jr-a)  ]  }  ; 


1  le 

1(n  -t  )  +  rP-^- 

2  ®  °°  1+a 


in 


where  e m  is  the  rotation  at  infinity,  N  is  the  complex  conjugate  of  N,  and 

T  ,  N  ,  R,  a  and  i  are  shown  in  Fig.  5.  To  obtain  the  solution  for  a 

constant  crack  surface  pressure  a,  which  is  the  loading  shown  in  Fig.  1(b) 
and  expressed  by  Eqs.  (1)  and  (2),  consider  the  stress  distribution  for  a 
cylindrical  inclusion  embedded  in  an  infinite  matrix  under  uniaxial  tension 
o 0  (Fig.  6).  The  radial  stress  along  the  fiber/matrix  interface  may  be 

expressed  as  [32]: 


°o 

<7rr  -  y(l-cos20) 


2a  o { A(^) 2  +  B [ 3 (~) 4 -4(^) 2 ] cos 20 ) ; 


(20) 
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where 


(l-2i/  )u  -  (1-2./  W 
.  _ t  m _ m  r 

4[(l-2 

and 


M  ‘M.C 
*m  't 


4[/i  +(3-4i/  )/if] 
m  m  ^fJ 


Superimposing  to  this  expression  a  uniaxial  loading  aQ  in  the  x  direction, 

with  r-R,  the  radial  stress  along  the  fiber/matrix  interface  becomes 
constant  (independent  of  6): 


arr  “  a°(1'4A)-  (21) 

Thus,  to  obtain  a  stress  cr-constant  at  the  interface,  we  must  have: 

a0  -  a/(l-4A).  (22) 

The  strain  energy  release  rate  at  the  crack  tip  for  the  thermal  loading 
is  then  obtained  by  using  Eq.  (19)  with  T^  ”  anci  Defining: 


GE 

2(1+1 /  ) a 2 R 
m 


(23) 


the  critical  fiber  size  for  the  suppression  of  interface  cracking  can  be 
obtained  by  equating  G  to  Gc ,  the  total  critical  strain  energy  release  rate 

of  the  interface.  Again,  noting  that  the  crack  length  is  not  known  a 

priori,  we  obtain  the  minimum  critical  fiber  size  R™ln  when  G'  is  maximum, 


G'  . 
max 


Thus : 
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(24) 


,mm 


G  E 
c  m 


2(l+i/  )ct2G  '  • 

m  max 


As  an  example,  we  assume  the  material  properties  given  in  the  previous 
section.  Figure  7  shows  the  variation  of  G'  with  the  crack  angle  a.  Again, 
it  is  seen  that  the  curve  reaches  a  maximum  of  G'  =*  0.53  at  a  =*  65°.  With 

these  values  and  AT  -  +1,000°C,  Eq.  (24)  gives  R^n  =.  H.7^m  which  is  of  the 

same  order  as  the  actual  fiber  size.  To  assess  the  effect  of  the  elastic 
properties,  the  ratio  of  Young's  moduli,  Em/E^,  is  varied  from  0.2  to  5. 

Unlike  the  case  of  radial  cracking,  as  Fig.  8  shows,  for  interface  cracking 
the  normalized  strain  energy  release  rate  G'  is  heavily  dependent  on  the 
elastic  properties  of  the  constituents. 

(b)  Interaction  of  two  fibers 

This  model  is  designed  to  simulate  situations  similar  to  those  seen  in 
Photo  2.  In  this  case  it  is  assumed  that  two  fibers  are  embedded  in  an 
infinite  matrix  and  subjected  to  mechanical  or  thermal  loads.  Again  we 
considered  two  crack  configurations:  (a)  a  radial  matrix  crack  and  (b)  an 
interface  crack  as  shown  in  Figs.  9(a)  and  9(b),  respectively.  Due  to  the 
complexity  of  the  geometry,  here  numerical  and  approximate  methods  have  been 
used  to  develop  the  solutions.  For  the  uncracked  geometry,  the  stress 
fields  in  the  matrix  and  at  the  fiber/matrix  interface  are  determined  using 
the  finite  element  technique.  We  have  developed  a  finite  element  computer 
procedure  using  ANSYS  [33]  which  gives  the  stress  and  displacement  fields 
for  in-plane  mechanical  and  thermal  loadings.  Figures  10  to  15  show  the 
typical  stress  distributions  of  the  two-fiber  model  due  to  thermal  loading 
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and  Figs.  16  to  21  show  the  stress  distributions  due  to  mechanical  loading. 
In  either  case  effects  of  various  Young's  modulus  ratios  and  fiber  spacings 
were  investigated.  It  is  found  that  the  stress  distributions  are  strongly 
influenced  by  fiber  interaction. 


(a)  Radial  Crack  in  Matrix  (b)  Interface  Crack 

FIGURE  9. 

The  crack  configurations  shown  in  Figs.  9(a)  and  9(b)  are  handled  as 
follows:  For  the  radial  crack,  we  developed  an  approximate  solution  in 
terras  of  a  singular  integral  equation.  The  approximation  is  due  to  the  fact 
that  the  kernel  used  in  the  singular  integral  equation  is  taken  from  the 
single  fiber  model,  while  the  crack  surfaces  are  loaded  with  the  negative  of 
the  stresses  obtained  above.  This  means  that  Eqs .  (5)  to  (18)  are  still 
valid  with  the  understanding  that  p(x)  in  Eq.  (5a)  is  now  replaced  by  the 
stress  obtained  for  the  uncracked  two-fiber  geometry.  The  numerical  results 
are  computed  for  varying  material  properties  and  fiber  separation  (25R)  and 
are  displayed  in  Figs.  22  to  24,  The  results  indicate  that  (a)  material 
properties  have  very  little  effect  on  the  normalized  fracture  energy,  g(«), 
as  shown  in  Fig.  22;  (b)  the  separation  plays  a  significant  role  on  whether 
crack  arrest  is  possible  (e.g..  Fig.  23  indicates  that  when  5  is  large, 
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i.e.,  the  fibers  are  far  from  each  other,  the  normalized  strain  energy 
release  rate  passes  through  a  maximum  making  crack  arrest  possible;  on  the 
other  hand  when  6  is  small,  i.e.,  when  the  fibers  are  close  to  each  other, 
there  is  no  possibility  of  radial  crack  arrest) ;  and  (c)  the  variation  of 
the  normalized  strain  energy  release  rate  g(e)  is  drastically  different  for 
the  single-fiber  and  the  two-fiber  geometries. 

For  the  interface  crack  problem,  again  an  approximate  solution  based  on 
Toya's  single-fiber  model  [31]  is  developed.  First  it  may  be  noted  that 
Toya's  solution  [31]  (Fig.  25(a))  can  be  considered  as  the  superposition  of 
an  uncracked  geometry  (Fig.  25(b))  and  a  crack  loaded  with  the  negative 
stresses  obtained  in  the  first  part  (Fig.  25(c)).  Then  the  strain  energy 
release  rate  can  be  computed  by  considering  only  the  perturbation  problem 
shown  in  Fig.  25(c).  Next  the  thermal-  or  mechanical -loading  induced 
tangential  and  radial  stresses  at  the  fiber/matrix  interface,  as  shown  in 
Figs.  10  to  21,  can  be  approximated  in  the  least  square  sense  to  yield  the 
crack  surface  tractions  prescribed  in  Fig.  5(c).  Thus  in  Toya's  solution, 
Eq.  (19),  the  forces  N  and  T  and  the  angle  tj>  are  replaced  by  N  ,  T  and 

<i>  ,  respectively.  The  numerical  results  are  displayed  in  Figs.  26  to  30 
eq 

tor  thermal  and  mechanical  loading,  where  the  normalized  strain  energy 
release  rate  is  plotted  against  the  interfacial  crack  angle  a.  The  material 
properties  and  the  crack  spacing  are  varied  and  the  results  are  compared 
with  those  obtained  for  the  single- fiber  case.  It  is  found  that  the 
properties  of  fibers  and  matrix  and  the  spacing  of  fibers  may  affect  the 
interfacial  crack  propagation  considerably.  Furthermore,  contrary  to  the 
radial  crack  case,  the  variation  of  G'  is  similar  qualitatively  for  both  the 
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single-fiber  and  the  two-fiber  cases,  making  interfacial  crack  arrest 
possible  (see  Figs.  27  to  30). 

(c)  The  ring  model 

This  model  is  designed  to  simulate  situations  similar  to  those  seen  in 
Photo  3.  The  model  consists  of  a  single  fiber  surrounded  by  a  circular 
array  of  fibers.  This  array  of  fibers  is  assumed  to  form  a  ring  of  material 
with  properties  obtained  by  considering  it  a  composite.  The  model  is  shown 
in  Fig.  31.  The  stress  distribution  in  the  uncracked  composite  due  to 
mechanical  and  thermal  loadings  can  be  obtained  in  closed  form  using 
Mitchell's  solution  [34]  for  the  four  regions  shown  in  Fig.  32. 


0  ©  '  ;'"0 

G/©0 1\ 

'  {  e-P, 


The  Ring  Model 


FIGURE  31. 


Assuming  plane  strain. 
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C  ,  -  €  ,  -  «  3  -  €  4~0; 

ZZ1  ZZ2  ZZJ  ZZ* 


'zz'  ■  E^zz'  •  +  *  °f4T; 


•5  ”  tr- [o  ,  -  v  (a  ~  +  a  AT; 

:z2  E  1  zz2  mv  rr2  86*  m 


m 


(25) 


(26) 

(27) 


18 


£  Z2  3 

E3  ^ azz3  1/3  (°rr3 

+  ^^^3)  ]  +  or3AT; 

(28) 

£zz4 

“E~[a2z4'  u*3  (arr 

4  +  +  a4  AT. 

(29) 

Equations 

(26) - (29)  can  also 

be  written  as : 

a  i 
zz1 

"  uf (orrx  +  a68 x) 

-  E^q^AT ; 

(30) 

a  2 
zz1 

"  "m^rr2  +  a6S^  ' 

•  E  a  AT ; 
m  m 

(31) 

a  , 
zz  3 

-  ^3(^rr3  +  °es 3)  ■ 

■  E3q3AT; 

(32) 

O  4 

zz4 

”  1/43  ^arr4+  a^4') 

■e43^at. 

(33) 

In  these  equations  Em,  Ef,  <*£  refer  to  the  Young's  moduli, 

Poisson's  ratios  and  thermal  expansion  coefficients  of  the  matrix  and 
fibers,  respectively,  while  u 3>  E43,i/43  are  the  composite  elastic 

properties  for  regions  3  and  4,  which  are  described  below  with  the  geometry 
shown  in  Fig.  33. 

Revlon  3  contains  a  group  of  fibers  randomly  distributed  in  a 
concentric  annular  domain  with  inner  radius  b  and  outer  radius  c.  The 

.  i’:r. :  modulus  E,  ,  r-'isson's  ratio  u,  and  thermal  expansion  coefficient  a, 

■j  J 

of  the  composite  are  dependent  on  n,  the  number  of  the  fibers  in  this  domain 

■  nu  tr. '!  “eometrical  parameters.  As  mown  below,  thev  can  be  determined 
virtue  of  the  rule  of  mixtures: 


E  -  E  -[n^-  +  E  c  — '  —  -  -  • 
ta  Lfc‘-b‘  +  m  c*-b2  ’ 


na :  .  c:  -  b  - - n  a  * 

J  fc*-b*  ra  C'-b2 


(24) 

(35) 
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(36) 


na2  c2-b2-na2 

°3  =  afc2 -b2  +  °m  c 2 ■ b 2  • 


Region  4  is  an  infinitely  extended  domain  with  a  hole  of  radius  c.  The 
material  in  this  region  is  transversely  isotropic  and  the  Young's  modulus, 
Poisson's  ratio  and  thermal  expansion  coefficient  of  the  material  are 
determined  by  the  following  formulas  (see  Fig.  33): 


V_  V  ’ 

-f+-E 

E*  E 
f  m 


where  E  is  the  Young's  modulus  in  Xj-  x2  plane  (see  Fig.  33),  is  the 
volume  fraction  of  fibers  defined  as: 


Vf  “ 


'n+2)ffa2 


and  V  -  1  -  V-; 

m  f  ’ 


with  n  the  number  of  fibers  inside  the  ring. 


"  -  Vf  +  (1-Vf^m  - 


4c2  -  (n+3~)7ra2 


Ac2  "f  + 


fn+3)*-a2 
— 7"n; — v  : 
4c2  m 


„  .  4c2  -  ( n+3 1  -rn2  v ..  ■  ^  , .. .» 

a  -  +  (1-Vf)am  -  4cz  a£  +  Ln2  am , 


fn+])7rq; 

4c2  am’ 


Uc^  -  ( n  +•  t  n  ^  ^ 

4c‘  4 c  *■  ’ 


where  v  and  a  are  the  Poisson's  ratio  and  the  thermal  expansion  coefficient 

in  the  plane  and  EtJ  is  th.e  Young's  modulus  in  the  x3  direction. 


Furthermore 
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=  a; 


(43) 


where  ^43  and  <*4  are  properties  in  x3  direction. 

Influence  of  the  surrounding  fibers  in  the  ring  model  is  considered  in 
the  average  sense.  The  influence  of  the  surrounding  fibers  can  be  adjusted 
either  by  changing  the  value  of  b,  i.e.,  changing  the  closeness  of  the  fiber 
cluster  or  by  changing  the  number  of  the  fibers  n  in  the  ring.  When  the 
number  of  the  fibers  in  the  ring  is  changed,  the  volume  fraction  of  the 
fibers  in  the  ring  is  changed;  and  thus  the  Young's  modulus  E3,  Poisson's 

ratio  i/3  and  the  thermal  expansion  coefficient  q3  will  all  change  and  the 

influence  of  the  ring  (fiber  cluster)  will  be  different. 

The  continuity  conditions  at  the  interfaces  can  be  written  as: 


U1 

“  u2  i 

vx 

-  v2, 

fr Bl  ~  rt9 2’ 

‘"rr1 

-  arr2’>  (r  “  a)  (44) 

u2 

-  u3, 

V2 

-  V3. 

rr 92  "  rr03’ 

‘"rr2 

“  arr3<  <r  “  b)  (45) 

U3 

-  U4, 

v3 

-  v4, 

rr93  ~  Tr 9*' 

CTrr3 

-  crrr4  .  (r  -  c)  (46) 

Solving 

the 

equations , 

we  obtain  the 

expressions  for  stresses  in  closed 

. ,  l :  J  v  H 1 

or. 

th.e 

•.press ion  which,  ar 

2  ver* 
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typical  distribution  of  agg  is  given  in  Fig.  34.  The  crack  problem  for  this 

L  j  :  icons oed  later  in  SvC .  5. 

4 .  Experimental  Work 

The  experimental  program  is  designed  to  accomplish  the  following:  (a) 
to  observe  and  document  the  raicrccracking  patterns  of  a  ceramic  -  fiber/ 
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ceramic-matrix  composite  under  high  temperature,  (b)  to  obtain  the  variation 
of  the  fracture  toughness  of  the  matrix  material  with  temperature,  (c)  to 
evaluate  the  effects  of  such  factors  as  fiber  size  and  distribution,  thermal 
expansion  coefficients,  and  temperature  rise  on  matrix  microcracking  of 
ceramic  composites,  and  (d)  to  determine  the  shear  strength  at  the 
fiber/matrix  interface. 

(a)  The  experimental  set-up 

The  experiments  consist  of  micro- indenting  monolithic  and  composite 
specimens  with  a  microhardness  indenter  under  various  temperatures  and 
recording  the  resulting  microcracking  patterns  using  a  microscope.  For  this 
purpose  an  ATS  series  3320  split- tube  laboratory  furnace  (Applied  Test 
Systems,  Inc.,  Butler,  PA)  which  can  be  heated  up  to  3000°F  in  30  minutes 
and  a  Nikon  UM-2  universal  measuring  microscope  are  used.  The  indentation 
load  is  recorded  by  a  10-lb  super-mini  load  cell  (Interface,  Inc., 
Scottsdale,  AZ)  and  the  temperature  is  read  by  a  type-B  platinum/platinum- 
30%  rhodium  thermocouple.  A  sketch  of  the  experimental  set-up  is  shown  in 
Fig.  35.  Figure  36  shows  the  series  3320  split- tube  laboratory  furnace. 

(b)  Effect  of  temperature  on  toughness  of  matrix  material 

The  toughness  of  the  matrix  material  is  obtained  by  indenting 
monolithic  SiC  coated  specimens  described  in  Sec.  3.  The  specimens  were 
first  mounted  in  a  bakelite  mold,  polished  with  paste  of  extra  fine  alumina 
and  diamond  powders  (as  small  as  half  a  micron)  and  indented  with  a  diamond 
Vickers  indenter  (Wilson  Instruments,  Inc.,  Bridgeport,  CT) .  They  were  then 
examined  under  optical  and  scanning  electron  microscopes  to  document  pre¬ 
existing  microcracks  and  only  a  few  of  them  were  observed.  To  ascertain  the 
eflect  of  temperature  on  the  toughness  of  the  material,  the  testing 
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temperature  was  varied  from  room  (25°C)  to  800°C.  Thus  micro-indentation 
tests  were  carried  out  at  25 °C,  250°C,  600°C  and  800°C.  At  each  temperature 
level  the  indentation  load  P  was  recorded  by  the  super-mini  load  cell  and 
the  half  size  of  the  impression  D,  as  shown  in  Fig.  37,  was  measured  by  the 
Nikon  UM-2  microscope.  From  this  information  one  may  calculate  the 
toughness  of  the  material  using  a  typical  formula  as  given  below: 


Kt  - 
Ic 


*3/2tanV>V/2 


), 


(47) 


where  P  is  the  applied  load,  D  the  half  crack  length  and  if>  -  68°  is  the 
indenter  angle.  Details  of  the  indentation  technique  can  be  found  in  Lawn 
et  al.  [19]  or  Evans  [20].  The  results  for  the  fracture  toughness  of  SiC 
at  various  temperatures  are  displayed  in  Fig.  38.  One  may  note  that,  as 
expected,  there  is  some  scattering  in  the  results.  However,  more 
importantly  it  seems  that  the  f-acture  toughness  of  SiC  decreases 
significantly  with  increasing  temperature.  For  example  KIc  decreases  from 

approximately  4.5  ksi/fn  at  room  temperature  to  3.5  ksi/In  at  800“C.  The 
fracture  toughness  of  the  matrix  was  also  measured  from  the  composite 
specimens.  To  preclude  any  effect  of  fibers  on  the  results,  the  indentation 
was  done  in  sites  with  no  neighboring  fibers.  The  same  results  were 
obtained,  as  will  be  seen  later  from  the  results  presented  in  Figs.  40-43. 
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electron  microscopes  to  document  the  pre-cracking  (existing  microcracks)  and 
the  fiber  distribution  patterns  used  in  the  analytical  model.  Again  only  a 
few  existing  microcracks  were  observed.  The  specimens  were  then  heated  in 
the  ATS  furnace  up  to  3000°F.  The  heated  specimens  were  later  examined 
under  optical  and  scanning  electron  microscopes.  Again  no  further 
microcracking  was  observed  due  to  temperature  rise  alone.  This  confirms  the 
claim  of  the  manufacturer  and  our  preliminary  calculations,  using  the 
single-fiber  model,  Sec.  3(a),  that  for  this  material  microcracks  do  not 
occur  under  thermal  loading  alone.  The  combination  of  thermal  and 
mechanical  loadings  is  necessary  to  generate  microcracks. 

Next  the  l/2"xl/2"xl/8n  Nicalcn/SiC  composite  specimens  with  fibers 
facing  up  were  also  polished  as  before  and  in  the  matrix  of  the  composite, 
microcracks  were  generated  again  using  the  indentation  technique  at  various 
temperatures  ranging  from  room  to  800°C.  To  ascertain  the  effects  of 
fibers,  the  specimens  were  indented  at  locations  of  varying  fiber  density 
and  the  fracture  toughness  was  calculated  using  Eqn.  (47),  as  if  the  fibers 
did  not  exist.  Since  the  fracture  toughness  described  above  contains  the 
effect  of  fibers,  from  here  on  it  will  be  referred  to  as  the  "apparent 
fracture  toughness"  of  the  matrix  material  (SiC)  .  The  effect  of  fibers  is 
introduced  through  the  concept  of  local  \rolume  fraction:  A  cell  of  fixed 
pronortions,  Fig.  39,  (sav,  b/a-3,  where  2a  is  the  crack  length  and  2b  is 
the  total  length  of  the  cell)  is  selected  with  the  indentation  imprint  at 
the  center  and  several  fibers  scattered  within  the  composite  cell.  The 
local  volume  fraction  is  defined  as  the  ratio  of  the  total  cross- 

sectional  areas  of  fibers  to  the  total  area  of  the  composite  cell.  As  a 
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consequence,  at  each  indentation  site,  one  obtains  a  different  local  V^. 

The  experimental  results  obtained  at  room  (25°C),  250°C,  600°C  and  800°C  are 
shown  in  Figs.  40  to  43.  In  each  figure  the  "apparent  fracture  toughness" 
is  plotted  against  the  local  volume  fraction  of  fibers,  V^,  at  a  given 

temperature.  It  is  observed  that  the  "apparent  fracture  toughness" 
decreases  with  increasing  local  volume  of  fibers  and  also  with  increasing 
temperature.  This  means  that  it  is  easier  to  generate  microcracks  at  a 
location  where  the  density  of  fibers  is  higher  and  when  the  temperature  is 
increased.  The  result  that  the  "apparent  fracture  toughness"  of  the  matrix 
is  decreasing  with  increasing  local  volume  fraction  must  be  expected,  since 
the  composite  under  consideration  is  a  so-called  "reversed"  composite, 
meaning  that  the  matrix  is  stiffer  than  the  fibers.  Here  the  Young’s 

modulus  of  the  matrix  material  (SiC)  is  Em~55xl0^  psi  and  the  Young's 

modulus  of  fibers  (Nicalon)  is  E^-26xl0^psi .  These  results  are  summarized 

in  Fig.  44,  where  the  experimental  data  are  presented  with  fitted  straight 
lines . 


(d)  Determination  of  interfacial  strength  and  interfacial  stresses 
The  interfacial  strength  at  debonding  initiation  is  determined  by  a 
combination  of  experimental  measurements  and  finite  element  analysis.  First 
the  load  at  debonding  initiation  is  measured.  For  this  purpose  some 
specimens  were  micro  -  indented  at  room  temperature  using  the  Wilson  Tukon  300 
microhardness  tester  whose  attached  microscope  enabled  us  to  indent  at  the 
desired  site  and  to  observe  the  ensuing  microcracks.  For  instance,  two 
perpendicular  matrix  cracks  generated  by  loading  a  Vickers  indenter  at  600 
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grams  is  shown  in  Photo  4.  Photo  5  shows  a  partial  interfacial  crack 
produced  at  the  interface  of  a  fiber  which  is  surrounded  by  a  circular  array 
of  fibers.  The  Vickers  indenter  was  loaded  at  50  grams  with  an  eccentricity 
of  approximately  three  quarters  of  the  fiber  radius.  Photo  6  depicts 
complete  debonding  of  two  isolated  fibers,  each  was  loaded  eccentrically 
with  a  Vickers  indenter  at  50  grams.  The  eccentricity  in  either  fiber  is 
about  one  half  of  the  fiber  radius.  From  these  observations  we  can  conclude 
that  the  micro- indentation  fracture  behavior  of  ceramic  composites  is 
affected  not  only  by  the  magnitude  and  location  of  the  indenter  but  also  by 
the  distribution  of  the  fibers,  i.e.  fiber  interaction. 

To  simplify  the  calculation,  the  indentation  is  performed  at  isolated 
fiber  locations  when  determining  the  interfacial  shear  strength.  On  can 
then  neglect  the  far  away  fibers  and  model  the  problem  as  a  fiber  surrounded 
by  an  infinite  matrix.  In  general,  however,  the  model  used  is  a  single 
fiber  surrounded  by  an  annulus  of  matrix  material  with  the  rest  of  fibers 
and  matrix  being  modeled  as  a  transversely  isotropic  composite  as  shown  in 
Fig.  45.  A  similar  model  was  developed  in  [35].  A  number  of  debonding 
experiments  were  conducted  and  in  each  case  the  critical  load  P  ,  which  is 

the  load  at  debonding  initiation,  and  the  fiber  size  were  recorded.  During 
the  debonding  test  the  indenter  is  presse  at  the  center  of  the  fiber  and 
the  lead  is  applied  incrementally  until  complete  debonding  occurs .  Through 
a  regression  analysis,  it  can  easily  be  shown  that  the  critical  loads  P 

are  proportional  to  the  square  of  fiber  radii  (see  Fig.  46).  This  finding 
indicates  that  the  results  are  consistent  for  all  the  fibers  with 
different  sizes. 
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To  relate  the  critical  debonding  load  to  the  strength  of  the  fiber- 
matrix  interface,  one  needs  the  distribution  of  interfacial  stresses.  For 
this  purpose  a  finite  element  model  has  been  developed.  In  the  finite 
element  model  the  loading  is  specified  as  a  concentrated  force  acting  at  the 
center  of  the  fiber  (Fig.  44(b)).  In  the  analysis  ANSYS  [33]  finite  element 
code  is  used.  Since  the  geometry  and  loading  are  axisymmetric ,  the 
axisymmetric  element  is  employed.  To  check  the  accuracy  of  the  solution, 
the  fiber  and  the  matrix  are  assigned  the  same  mechanical  properties  and  the 
results  obtained  are  compared  with  the  closed  form  Boussinesq's  solution 
[34];  The  results  agree  extremely  well.  The  calculations  are  performed 
in  two  steps,  i.e.  first  a  course  mesh  is  used  and  later  a  finer  mesh  is 
employed,  A  typical  mesh  used  in  the  computation  is  shown  in  Fig.  46ct.  The 
normalized  shear  and  radial  stresses  at  the  fiber-matrix  interface  are  shown 
in  Fig.  47.  The  stresses  due  to  temperature  change  were  also  considered  and 
are  displayed  in  Fig. 48.  The  debonding  shear  strength  is  the  maximum  shear 
strength  when  the  load  is  Pcr-  Based  on  the  analysis  given  above,  the  shear 

strength  can  be  calculated  by  the  following  formula: 


P  /R-’ 
cr 


K _ 

1 

+  E  q  AT 
m  m 


(48) 


where  K  is  a  constant  (in  our  example  K  =  0.1S),  Pc>_  the  critical  load 

(i.e.,  the  load  at  the  onset  of  debonding),  E  Young's  modulus  of  the 

ra 

matrix,  a  thermal  expansion  coefficient  of  the  matrix,  R  radius  of  the 
m 

fiber  being  tested  and  AT  change  in  temperature. 
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5 .  Theoretical  Prediction  and  Comparison  with  Experimental  Results 

To  explain  and  predict  the  observed  experimental  results  given  in  the 
previous  section,  a  theoretical  model  has  been  developed.  In  a  composite 
material,  the  distribution  of  fibers  may  assume  many  patterns  such  as  an 
isolated  fiber,  two  interacting  fibers,  a  ring  of  fibers,  or  random.  Here 
we  assume  that  the  fibers  form  a  ring  of  composite  material  with  elastic 
properties  calculated  based  on  a  volume  fraction  equivalent  to  the  volume 
fraction  in  the  cell  described  above.  The  geometry  of  the  model  is  shown  in 
Fig.  49.  Since  the  properties  of  the  constituents  are  fixed,  from 
dimensional  considerations,  it  can  be  shown  that  for  the  ring  model  the 
stress  intensity  factor  for  a  crack  of  length  2a  in  a  uniformly  loaded 
composite  can  be  expressed  as: 

Kj  -  ajna.  f(Vf,  R/a,  AT),  (49) 


where  a  is  the  crack  surface  traction,  V-.  the  local  fiber  volume  fraction, 
o  f 

R  the  fiber  radius,  and  AT  the  temperature  rise.  When  oq  reaches  the 

critical  value  (<7q)c,  will  become  the  toughness  of  the  matrix, 

(K,  )  .  .  Thus  the  "apparent  fracture  tounhness",  which  is  defined  as: 

ic  matrix 


(Kt  )  -  (a  )  /ira , 

Ic  apparent  o  c 


can  be  expressed  in  terms  of  and  R/a  as: 
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(52) 


f(V£,  R/a,  AT)  -  7^7l, 

o 

is  the  normalized  stress  intensity  factor. 

The  stress  intensity  factor  Kj.  is  determined  using  the  finite  element 

program  ABAQUS  [36].  Thus  for  a  given  R/a  and  AT,  a  (Kic)apparent  versus  Vf 

curve  can  be  plotted.  Numerical  results  for  R/a  -  0.95  and  T  -  25°C,  250°C, 
600°C,  and  800°C  have  been  computed.  The  finite  element  mesh  used  in  the 
calculations  is  shown  in  Fig.  50.  The  theoretical  results  obtained  for  each 
temperature  level  are  plotted  against  the  experimental  results  and  are 
depicted  in  Figs.  51  to  54.  It  is  seen  that  the  agreement  between  the  two 
sets  of  results  is  good.  The  model  also  predicts  a  decreasing  "apparent 
fracture  toughness"  with  local  volume  fraction  of  fibers.  For  low  volume 
fractions  and  high  temperatures,  the  quantitative  prediction  is  excellent. 
However,  for  large  volume  fractions  and  low  temperatures,  there  is  some 
deviation  between  the  experimental  and  the  predicted  results.  This  is 
expected  since  our  model  does  not  take  into  account  the  specific 
arrangement  of  fibers  around  the  crack  tip  and  the  possible  mode  II  effects, 
which  mav  influence  the  "apparent  fracture  toughness"  significantly. 

6 .  References 

1  ’  A.  G.  Evans  and  K.  T.  Faber.  "Crack-Growth  Resistance  of 
Microcracking  in  Brittle  Materials,"  Journal  of  the  American 
Ceramic  Society.  Vol.  67,  pp .  255-260,  1984. 

[2]  R.  W.  Rice,  "  Fractographic  Identification  of  Strength  Controlling 
Flaws  and  Microcracking,  :  in  Fracture  Mechanics  of  Ceramics.  Vol. 
1,  ed.  R.  C.  Bradt,  D.  P.  H.  Hasselman  and  F.  F.  Lange,  Plenum 
Press,  New  York,  1974,  pp .  323-345. 


29 


[3]  R.  W.  Rice,  S.  W.  Frieraan,  R.  C.  Pohanka,  J.  J.  Mecholsky,  and  C. 
C.  Wu,  "Microstructural  Dependence  of  Fracture  Mechanics 
Parameters  in  Ceramics,"  in  Fracture  Mechanics  of  Ceramics.  Vol. 
4,  ed.  by  R.  C.  Bradt,  D.  P.  H.  Hasselman  and  F.  F.  Lange,  Plenum 
Press,  New  York,  1978,  pp .  840-876. 

[4]  R.  W.  Rice,  "Microstructural  Dependence  of  Mechanical  Behavior  of 
Ceramics,"  in  Treatise  on  Materials  Science  and  Technology,  Vol. 
2.  Properties  and  Microstructure,  ed.  by  R.  K.  McCrone,  Academic 
Press,  New  York,  1977,  pp .  199-381. 

[5]  R.  W.  Rice,  R.  C.  Pohanka  and  W.  J.  McDonough,  "  Effect  of 
Stresses  from  Thermal  Expansion  Anisotropy,  Phase  Transformations 
and  Second  Phases  on  the  Strength  of  Ceramics,"  Journal  of  the 
American  Ceramic  Society.  Vol.  63,  pp.  703-710,  1980. 

[6]  R.  W.  Rice,  "Processing  Induced  Sources  of  Mechanical  Failure  in 
Ceramics,"  in  Processing  of  Crystalline  Ceramics.  ed.  by  H. 
Palmour,  III,  R.  F.  Davis  and  T.  M.  Hare,  Plenum  Press,  New  York, 
pp.  303-319,  1978. 

[7]  F.  I.  Baratta,  "Stress  Intensity  Factor  Estimates  for  a 
Peripherically  Cracked  Spherical  Void  and  a  Hemispherical  Surface 
Pit,"  Journal  of  the  American  Ceramic  Society.  Vol.  61,  pp.  490- 
493,  1978. 

[8]  A.  G.  Evans,  D.  R.  Biswas  and  R.  M.  Fulrath,  "Some  Effects  of 

Cavities  on  the  Fracture  of  Ceramics:  I,  Cylindrical  Cavities," 
Journal  of  the  American  Ceramic  Society.  Vol.  62,  pp.  95-100, 
1979. 

[9]  A.  G.  Evans,  D.  R.  Biswas  and  R.  M.  Fulrath,  "Some  Effects  of 

Cavities  on  the  Fracture  of  Ceramics:  II,  Spherical  Cavities," 
Journal  of  the  American  Ceramic  Society.  Vol.  62,  pp .  101-106, 
1979. 

[110  D.  J.  Green,  "Stress  Intensity  Factors  Estimates  for  Annular 

Cracks  at  Spherical  Voids,"  Journal  of  the  American  Ceramic 

Society.  Vol.  63,  pp .  342-344,  1980. 

[11]  F.  I.  Baratta  and  A.  P.  Parker.  "Stress  Intensity  Factor  Estimates 

for  Various  Configurations  Involving  Cracked  Spherical  Voids,"  in 
Fr  v  -ur  2  ' '  ■  v  rni  £  ■  r  ■-  ;  ■-  s  .  Vol.  5,  ed.  by  R.  C.  Bradt,  A.  G. 

Evans,  D.  ?.  H.  Hasselman  and  F.  F.  Lange,  Plenum  Press,  New  York, 
1983,  pp.  543-567. 

[12]  R.  W.  Rice,  "Test-Microstructural  Dependence  of  Fracture  Energy 
Measurements  in  Ceramics,"  in  Fracture  Mechanics  Methods  for 
Ceramics .  ASTM  Symposium  on  Fracture  Mechanics  Methods  for 
Ceramics,  Rocks  and  Concrete,  Chicago,  1980. 


30 


[13]  R.  W.  Davidge  and  D.  J.  Green,"  The  Strength  of  Two-Phase  Ceramic- 
Glass  Material,"  Journal  of  Materials  Science.  Vol.  3,  pp .  629- 
634,  1968. 

[14]  F.  F.  Lange,  "Criteria  for  Crack  Extension  and  Arrest  in  Residual 
Localized  Stress  Field  Associated  with  Second  Phase  Particles,"  in 
Fracture  Mechanics  of  Ceramics.  Vol.  2,  ed.  by  R.  C.  Bradt,  D.  P. 
H.  Hasselman  and  F.  F.  Lange,  Plenum  Press,  1974 

[15]  A.  G.  Evans,  "The  Role  of  Inclusions  in  the  Fracture  of  Ceramic 
Materials,"  Journal  of  Materials  Science.  Vol.  9,  pp.  1145-1152, 
1974. 

[16]  D.  J.  Green,  "Stress - Induced  Microcracking  at  Second  Phase 
Inclusion,"  Journal  of  the  American  Ceramic  Society.  Vol.  64,  pp. 
138-141,  1981. 

[17]  Y.  M.  Ito  and  R.  B.  Nelson,  "Numerical  Modeling  of  Microcracking 
in  Two-Phase  Ceramics,"  in  Fracture  Mechanics  of  Ceramics.  Vol.  5, 
ed.  by  R.  C.  Bradt,  A.  G.  Evans,  D.  P.  H.  Hasselman  and  F.  F. 
Lange,  Plenum  Press,  New  York,  1983,  pp.  479-493. 

[18]  Y.  M.  Ito,  M.  Rosenblatt,  L.  Y.  Cheng,  F.  F.  Lange,  and  A.  G. 
Evans,  "Cracking  in  Particulate  Composites  due  to  Thermomechanical 
Stress,"  International  Journal  of  Fracture.  Vol.  27,  pp.  483-491, 
1981. 

[19]  B.  R.  Lawn  and  T.  R.  Wilshaw,  "Review  -  Indentation  Fracture: 
Principle  and  Applications,"  Journal  of  Materials  Science.  Vol. 

10,  pp.  1049-1081,  1975. 

[20]  A.  G.  Evans,  "Fracture  Toughness:  The  Role  of  Indentation 
Technique,"  Fracture  Mechanics  Applied  to  Brittle  Materials.  ASTM 
STP  678,  S.  W.  Freiman,  Ed.,  American  Society  for  Testing  and 
Materials,  pp.  112-135,  1979. 

[21 1  F.  Erdor.an  and  G.  D.  Gupta.  "The  Inclusion  Problem  with  a  Crack 
Crossing  the  Boundary",  International  Journal  of  Fracture.  Vol. 

11,  pp.  13-27,  1975. 

[22]  F.  Erdogan.  "Mixed  Boundary  Value  Problems,"  in  Mechanics  Today, 
Vo],  4.  ed.  bv  S.  Nemat-Nujser ,  Pergamon  Press,  Oxford,  1978,  pp . 
1-2-;. 

[23]  D.  B.  Marshall  and  A.  G.  Evans,  "Failure  Mechanisms  in  Ceramic- 
Fiber/Ceramic  -Matrix  Composites",  Journal  of  the  American  Ceramic 
Soc iety .  Vol.  68,  pp.  225-231,  1985. 

[24]  F.  Erdogan,  G.  D.  Gupta  and  T.  S.  Cook,  "Numerical  Solution  of 
Singular  Integral  Equations,"  Mechanics  of  Fracture.  Vol.  1: 
Methods  of  Analysis  and  Solutions  of  Crack  Problems,  ed.  by  G .  C. 
Sih,  Noordhoff,  Leyden,  The  Netherlands,  1973,  pp .  368-425. 


31 


[25]  N.  I.  Iaokimidis  and  P.  S.  Theocaris,  "On  the  Solution  of 
Collocation  Points  for  the  Numerical  Solution  of  Singular  Integral 
Equations  with  Generalized  Kernels  Appearing  in  Elasticity 
Problems,"  Computers  and  Structures.  Vol.  11,  pp.  289-295,  1980. 

[26]  A.  H.  England,  "A  Crack  Between  Dissimilar  Media,"  Journal  of 
Applied  Mechanics.  Vol.  32,  pp .  400-402,  1965. 

[27]  F.  Erdogan,  "Stress  Distribution  in  Bonded  Dissimilar  Materials 
with  Cracks,"  Journal  of  Applied  Mechanics.  Vol.  32,  pp.  403-410, 
1965. 

[28]  F.  Erdogan  and  G.  D.  Gupta,  "The  Stress  Analysis  of  Multi-Layered 
Composites  with  a  Flaw,"  International  Journal  of  Solids  and 
Structures .  Vol.  7,  pp .  39-61,  1971. 

[29]  M.  L.  Williams,  "The  Stress  around  a  Fault  or  Crack  in  Dissimilar 
Media,"  Bulletin  of  the  Seismological  Society  of  America.  Vol.  49, 
pp.  199-204,  1959. 

[30]  G.  P.  Cherepanov,  Mechanics  of  Brittle  Fracture.  McGraw-Hill,  New 
York,  1979. 

[31]  M.  Toya,  "A  Crack  Along  the  Interface  of  a  Circular  Inclusion 
Embedded  in  an  Infinite  Solid",  Journal  of  the  Mechanics  and 
Physics  of  Solids.  Vol.  22,  pp.  325-348,  1974. 

[32]  J.  N.  Goodier,  Journal  of  Applied  Mechanics.  Vol.  1,  p.  39, 

1933. 

[33]  ANSYS,  Swanson  Analysis  System,  Inc. 

[34]  R.  W.  Little,  Elasticity.  Prentice -Hall ,  Englewood  Cliffs,  NJ , 
1973. 

[35]  D.  H.  Grande,  J.  F.  Mandell  and  K.  C.  C.  Hong,  "Fibre -Matrix  Bond 

Strength  Studios  of  Glass,  Ceramic,  and  Metal  Matrix  Composites," 
Journal  of  Materials  Science.  Vol.  23,  pp .  311-328,  1988. 

[36]  ABAQUS ,  Hibbit,  Karlsson  &  Sorensen,  Inc.,  Providence,  RI . 


32 


J 


Fiber  Surrounded  by  fln 
Qrralf  cf-  fibers 


Photo  4.  hatux  traces 


Partial  Debonding  Along  I'M  k  4  Complete  Debonding  Along 
Fiber/Matrix  Interface  Fiber/Matrix  Interface 


38 


<jr 


1 


•o  «n  *'» 

Q' - n _ jOl _ Ol _ . _ 0_ 


matrix 


ation  of  G  with  the  interface  crack  angle  o(  for  various 
es  of  the  ratio  E^/E,. 


•))JLU|1UJJ  »m  o)  «jjidUJ!i!!K  0: 

- h, _ _ _ _ _ _ _ _ _ _ _ _ _ . _ _  _ _ 


>  i  ■  i  ho  ( '.»■ iiitnii »«  r 


■i’igurfi  15. 


OF  MATERIAL  PROPERTIES 

ITERFACIAL  CRACK(TWO  FIBER  MODEL) 


80  100  120  140 

CRACK  ANGLE  (DEGREES) 

O  Em/Ef  =  1  +  Em/Ef=2.2 

Figure  26. 


ITH  SINGLE  FIBER  MODEL 

INTERFACIAL  CRACK 


SINGLE  FIBER  MODEL 


;t  of  neighboring  fiber 

INTERFACIAL  CRACK  (TWO  FIBER  MODEL) 
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INTERFACIAL  CRACK 
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TWO  FIBER  MODEL  +  SINGLE  FIBER  MODEL 

Figure  30. 


a)  Calculation  of  volume  fraction  b)  Calculation  of  overal  volume 

of  fibers  in  region  3  fraction  in  region  4 


c  Coordinate  system  for  calculation 
cf  elastic  parameters 


tic  parameters 


fig  33 . Calculation  of  e' 
in  ring  model 
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VOLUME  FRACTION  Vf 


K1C  .vs.  Vf 

TEMPERATURE  250C 


72 


(0 


(  HLk  d(3r>o  ix 


VOLUME  FRACTION  Vf 


K1C  .vs.  Vf 

TEMPERATURE  600C 
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VOLUME  FRACTION  Vf 


VOLUME  FRACTION  Vf 


WITH  TEMPERATURE  EFFECT 
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VOLUME  FRACTION  Vf 


Figure  45. 


Debonding  load  .vs.  Fiber  Area 


Fig . 47  .Stress  distributions  on  the  interface 

along  the  fiber  axial  direction  under  debonding 
loading 
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Volume  Fraction  Vf 
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